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ABSTRACT 

In the recent paper of iHooper & Goodenoughl (1201 Ol) it was reported that 7-ray emis- 
sion from the Galactic Center region contains an excess compared to the contributions 
from the large-scale diffuse emission and known point sources. This excess was argued 
to be consistent with a signal from annihilation of Dark Matter with a power law den- 
sity profile. We reanalyze the Fermi data and find instead that it is consistent with the 
"standard model" of diffuse emission an d of known point sources . The m ain reason for 
the discrepancy with the interpretation of lHooper & Goodenoughl (|2010!) is different (as 
compared t o the previous works) sp e ctrum of the point source at the Galactic Center 
assumed by iHooper & Goodenoughl (l2010h . We discuss possible reasons for such an 
interpretation. 



1 INTRODUCTION 



The origin of the emission from the Galactic Center (GC) 
at keV-TeV energies has been extensively discussed in 
the literature over las t few years. In their recent paper, 
iHooper & Goodenouglj ilOld) claimed that the 7-ray emis- 
sion from the Ga lactic Center region , measured with the Fermi 
LAT instrument jAtwood et ai]|2009h cannot be described by a 
combination of spectra of known point sources, diffuse emis- 
sion from the Galactic plane and diffuse spherically symmet- 
ric component (changing on the scales much larger than 1°). 
An additional spherically symmetric component was suggested 
to be needed in the central several degrees. This component 
was then interpreted as a dark matter annihilation signal with 
the dark matter distribution having power law density profile 
p{r) (X r~", Q ~ 1.34. The observed excess is at energies be- 
tween ~ 600 MeV and ~ 6 GeV and the mass of the proposed 
DM particle was suggested to be in the GeV energy band. 

In this work we anal yze the Fermi data, used in 
iHooper & Goodenoughl JloToh . utilizing the data analysis tool, 
provided by the Fermi team. 



2 DATA 

For our analysis we consider 2 years of Fermi data collected be- 
tween August, 4th, 2008 and August 18th, 2010. The standard 
event selection for source analysis, resulting in the strongest 
background-rejection power {diffuse event class) was appliedQ 
In addition, photons coming from zenith angles larger than 105° 



were rejected to reduce the background from gamma rays pro- 
duced in the atmosphere of the Earth. 

The Fermi's point-spread function (PSF ) is non- 
Gaus s ian and strongly depends on energy jAbdo et afl 
I2OO9I : I Atwood et al. II2OO9'). In order to properly take it into 
account and better constrain the contributions from Galactic 
and Extragalactic diffuse backgrounds we analyze a 10° x 10° 
region around the Galactic Center. 



2.1 Model 

To describe emission in the 10° x 10° region we use the model 
containing two components - point sources and diffuse back- 
grounds. 

To model the contribution from the point sou rces we in- 
clude 19 sources from 11 months Fermi catalog jAbdo et aP 
l2010j) falli ng into the selected region plus 4 additional sources 
described in Chemvak ova et al] boioh . We fix the positions of 
the sources to coordinates given in the cata log. We model thei r 
spectra as power law (in agreement with lAbdoetaljr2010ah . 
Thus we have 46 free parameters (power law index and norm 
for each of the sources) to describe the point-source component 
of the model. 

To describe the diffuse component of emission, we use the 
models for the Galactic diffuse emission (gll_iem_v02 . fit) 
and isotropic (isotropic_iem_v002.txt) backgrounds 
that were developed by the LAT team and recommended for 
the high-level analysis j Abdo et ^ l2010blFI . These models 
describe contributions from galactic and extragalactic diffuse 



See e.g. http : //f ermi . gsf c . nasa.gov/ssc/data/ 
analysis/scitools 



|http : //fermi . gs f c .nasa . gov/ssc/data/analysis/] 
scitools /likelihood_tutorial . html 
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Figure 1. The map of significance of residuals for the region around the 
Galactic Center. 



backgrounds correspondingly. The number of free parameters 
for the diffuse background model is 2 (the norms for each of the 
backgrounds). The total number of free parameters in our model 
is thus 48. 

This model is s imilar to the one described in 
IChemvakova et aljfeoioh . 
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Figure 2. Spectrum of the point source at the GC reported in 
IChernvakova et all j2010l) (green points) together with the HGIO total 
spectrum from 1.25° (black points), excess (blue squares) and GC point 
sourc e flux from HGIO (red open circles). Continu ation of the HESS 
data ivan Eldik etai]l20O8l : lAharonian et alj|2004l) (blue points) data 
with a power law is shown with dashed black line. 



2.2 Analysis 

The data analysis was performed using the LAT Science Tools 
pack age with the P6-V 3 post-launch instrument response func- 
tion teando et alj2009l) . 

We find the best-fit values of all parameters of the model 
of Section|2T|(using qtlike likelihood fitting to ol) and deter- 
mine resulting log-likelihood ( iMattox et al.ll99q) of the model. 
Best fit values for the obtained fluxes agree withi n statistical 
uncert ainties w ith fluxes reported i n Ferm i Catalog jAbdo et alj 
l2010ah and in IChernvakova et al. I ( l2010h (e.g. for the central 
source we obtained the flux 5.68 x 10~* cts/cm^/s while the 
catalog gives (5.77 ± 0.3) x 10"* cts/cm^/s). 

We then freeze the values of the free parameters of our 
model and simulate spatial distribution of photons at energies 
above 1 GeV (using gtmodel tool). The significance of resid- 
uals, (Observation - Model)/ statistical error, is shown in Fig[T] 
We see the absence of structures in the central 2° region. The 
average value of residuals is about 10% in the 2° region around 
the GC, compatible with estimated systematic errors (10-20%) 
of Fermi LAT at 1 GeV0 

Thus we see that the adopted model (point sources plus 
galactic and extragalactic diffuse components) explains the 
emission from the GC region and no additional components is 
required. 



3 DISCUSSION 

We conclude that the si gnal within central 1°— 2°, c ontain- 
ing the "excess" found bv lHooper & Goodenouehlboid (HGIO 
hereafter), can be well described by our model : (point sources 
plus Galactic and extragalactic diffuse background compo- 
nents). The discrepancy is then due to a different interpretation 
of the data. 

^ See e.g. |http : //f ermi . gsf c .nasa ■ gov/ssc/data/ 
analysis/LAT_caveats .html^ 



The spectrum of the central point source (IFGL J1745.6- 
2900c, probably associated with the Galactic black hole Sgr 
A*) was taken in HGIO to be a featureless power-law start- 
ing from energies about 10 TeV (results of HESS measure- 
ments, blue points with error bars in Fig. (2] iAharonian et alj 
2004; van Eldik et al. 2008)) and continuing all the way down 
to ~ 1 GeV. The flux attributed in this way to the central 
point source is significantly weaker than in the previous works. 
For comparison, t he (PSF corrected) s pectru m of the GC point 
source reported in IChernvakova et alj j201ol) is shown in Fig.|2] 
in green points. Its spectral characteristics are fully consisten t 
with the results of 1 1 -months Fermi catalog lAbdo et al] i2010ah 
(~ 6 X 10~* cts/cm^/s above 1 GeV, compared to the ~ 
5 X 10~^ cts/cm^/s at the same energies in HGIO). The change 
of the slope of the source spectrum below ~ 100 GeV, as com- 
pared with the HESS data is explained by Chernvak ova et al] 
( l2O10h with the model of energy dependent diffusion of pro- 
tons in the few central parsecs around the GC. Alternatively, 
the spectrum can be explain ed with the model developed in 
lAharonian & Neronovl ( l2005 l). The low-energy (GeV) compo- 
nent of the spectra in this model is explained by synchrotron 
emission from accelerated electrons, while high-energy (TeV) 
one by inverse Compto n radiation of the same particles. Accord- 
ing to the analysis of lAbdo et al.l bOlOah : Ichemvakova et al] 
( l2010l) the central point source provides significant contribution 
to the flux in the 1 .25° central region. HGIO suggest, apparently, 
a different interpretation. They assume that there is no signifi- 
cant change in the spectrum of the central source at ~ 100 GeV 
and the spectrum observed by HESS at high energies continues 
to lower energies. Then, large fraction of the flux between the 
energies ~ 600 MeV and ~ 6 GeV has to be attributed to the 
"DM excess". One of the reasons in favor of such an interpreta- 
tion could be the feature in the total spectrum from the central 
region (rise between ~ 600 MeV and several GeV) discussed in 
HGIO. Such a feature would also be c onsistent with a possible 
contribution from millisecond pulsars (lAbazaiianll2oTob . that is 
also expected to have a maximum at ~ 2 — 3 GeV. 

To illustrate the nature of the spectral shape at these ener- 
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Figure 3. Left: the "inner" (5° ai'ound the Galactic plane) and "outer" regions. Right: Effects of the energy dependence of the effective area for the 
spectra of the "inner" and "outer" regions. 



gies we collected "front converted" (FRONT) photons from the 
region of the width 5° around the Galactic Plane (the ''inner" 
region) and from the ''outer" region as demonstrated on the left 
panel in Fig. [3] The count rate from each of these regions was 
divided by the constant effective area (3500 cm"^) to obtain the 
fluxQone sees that the total emission from both regions demon- 
strates the same spectral behavior as the excess of HGIO, sug- 
gesting that this spectral shape is not related to the physics of 
the several central degrees. This drop of flux at low energies is 
mainly due to the decreasing effective area of the satellite[f| If 
we properly take into account the dependence of the effective 
area on energy, we obtain the spectrum that "flattens" at small 
energies and exceeds by a significant factor the flux from the 
central point source (as it should) (compare red and magenta 
points on the right panel in Fig.[3)- 

Another reason for the decrease of the HGIO spectrum 
is the increase of Fermi LAT PSF at low (< 1 GeV) ener- 
gies[f| This means that if one collects photons from a relatively 
small region, such that a contribution from its boundary (with 
the PSF width) is comparable to the flux from the whole re- 
gion, the spectrum would artificially decline, due to increas- 
ing loss of photons at low energies. To disentangle properly 
what photons in the PSF region had originated from a localized 
source, and what are parts of the diffuse background, special 
modeling is needed. In the monoto nic spectrum of the GC, ob- 
tained bv lChernvakova et al.l ( 1201 Ol) both these effects (effective 
area and PSF) were taken into account as it was obtained from 
10° X 10° region, using the Fermi software. 

To further check the nature of the emission from the cen- 
tral several degrees, we took a fiducial model, that contained 
the same galactic and extragalactic diffuse components plus all 
the same point sources, but excluding the point source in the 
center. We then fit our data to this new model. Such a fit at- 



The effective area of Fermi LAT is strongly energy dependent. The 
number 3500 cm^, roughly corresponding to the effective area at r^ 
1 GeV, is used here as a quick expedient (see below). 
^ |http : //www-glast ■ slac ■ stanf or d ■ edu/ software/ 
IS/glast_lat_perf ormance -htm 

® For example, for normal incidence 95% of the photons at 1 GeV are 
contained within ~ 1.6° and in 2.8° at 500 MeV 



tempts to attribute as many photons as possible from the region 
around the GC to the emission of diffuse components. The pro- 
cedure leaves strong positive residuals within the central 1 — 2°. 
The spectrum of these residuals is consistent with the spectrum 
of the central point source of Chernvakova et al. (2010) (green 
points in Fig. [JJ. To demonstrate, that the spatial distribution 
of these residuals is fully consistent with the PSF of Fermi, we 
compare their radial distribution in various energy bins with the 
radial distribution around the Crab pulsar (as it was done e.g. 
in Neronov et al. (2010)). The pulsar win d nebu la, associated 
with the Crab has an angular size ~ 0.05° jHesteil2 008). Thus, 
for Fermi LAT Crab is a point source. The radial profile of resid- 
uals at all energies has the same shape as Crab, as Fig.|4]clearly 
demonstrates. As an additional check, we repeated the above 
test using only FRONT photons (as in this case the PSF is more 
narrow) and arrived to the same conclusion. 



The above analysis demonstrates that the emission around 
the GC in excess of diffuse components (galactic and extra- 
galactic) is fully consistent with being produced by the point 
source with the power-law spectrum, obtained in lAbdo et ahl 
( l2010ah :l Chernvakova et al ] (f20im . and no additional compo- 
nent is required. 



A different question however is whether such an addi- 
tional component may be ruled out. To this end we have added 
to our model of Sec l2. II an additional spherically symmetric 
component, whose intensity is distributed around the center as 
p (r) (where p{r-) oc as found in HGIO). We observe, 

that such a procedure does improve the fit (change in the log- 
likelihood is 25 with only one new parameter added). The re- 
sulting spectral component is shown in Fig.|5] Some of the pho- 
tons from the galactic diffuse background were attributed by 
the fit procedure to the new component, concentrated in several 
central degrees (within the Galactic Plane). This phenomenon 
is probably related to the complicated and highly non-uniform 
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Figure 4. Radial profile of residuals at different energies around the GC as compared to the radial profile of Crab emission (renormalized so that the 
total flux in each energy range coincide). In both cases only FRONT photons were used. 
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Figure 5. Spectrum of an additional spherically symmetric component, 
distributed around the GC as the HGIO excess. 

in the central region galactic diffuse backgrouncQ (cf. also the 
right panel of the Fig.|6}. 

^ See "Description and Caveats for the LAT Team Model of 
Diffuse Gamma-Ray Emission" by the Diffuse and Molecu- 
lar Clouds Science Working Group, Fermi LAT Collaboration, 
http : / / f erini . gsf c ■ nasa ■ gov/ssc/data/ access/1 at/ 
ring_f or_FSSC_f inal4 .pdfl 



We should also note that HGIO modeled diffuse back- 
ground differently. They considered contributions from the 
Galactic disk and spherically symmetric emission in the re- 
gion outside central 2° and then extrapolated the diffuse model 
into the innermost 1° — 2°, arguing that the contribution does 
not vary significantly in the range 2 ° — 10° off - center. The 
background model we used (see Abdo et al ] |2010j : lAbdoet alj 
l2010bl for the detailed description) is different from that of 
HGIO, especially in the central 1-2°, where the model flux is 
higher than the one extrapolated from larger galactic longitudes, 
as one can clearly see on the right panel of the Fig.|6] 



Having the above considerations in mind, we think that 
the spectrum of the central region, changing monotonously with 
the energy, is well described by purely astrophysical model of 
the central point source and therefore present data do not re- 
quire any additional physical ingredients, such as DM annihi- 
lation signal or additional contributions from millisecond pul- 
sars. However, to firmly rule out the emission from DM anni- 
hilation in the GC, more detailed model of the galactic diffuse 
background is required. Additionally, with the future data, bet- 
ter statistics will reduce the error bars on the data point around 
~ 100 GeV which will be helpful to better understand the cen- 
tral point source physics. 
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Figure 6. Left: 10°xlO° count map of best-fit model. Right: only contribution from galactic and extragalactic backgrounds is shown 
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